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Abstract 

OPERA has claimed the discovery of superluminal propagation of 
neutrinos. We analyze the consistency of this claim with previous 
tests of special relativity. We find that reconciling the OPERA 
measurement with information from SN 1987a and from neutrino 
oscillations requires stringent conditions. The superluminal limit 
velocity of neutrinos must be nearly flavor independent, must de- 
crease steeply in the low-energy domain, and its energy dependence 
must depart from a simple power law. We construct illustrative 
models that satisfy these conditions, by introducing Lorentz viola- 
tion in a sector with light sterile neutrinos. We point out that, quite 
generically, electroweak quantum corrections transfer the informa- 
tion of superluminal neutrino properties into Lorentz violations in 
the electron and muon sector, in apparent conflict with experimen- 
tal data. 



1 Introduction 



The OPERA collaboration has recently announced the measurement of the velocity of neu- 
trinos, as they travel from CERN to the Gran Sasso Laboratory (LNGS) covering a distance 
of about 730 km [lj. The CERN Neutrino beam to Gran Sasso (CNGS) consists of u^, with 
small contaminations of (2.1%) and of v e or v e (together less than 1%). The average 
neutrino energy is 17 GeV. OPERA has measured the neutrino velocity by taking the ra- 
tio between very accurate determinations of distance and time of flight. The distance is 
defined as the space separation between the focal point of the graphite target (where the 
proton beam extracted from the CERN SPS collides producing secondary charged mesons 
that eventually decay into neutrinos) and the origin of the OPERA detector reference frame. 
High-accuracy GPS measurements and optical triangulations led to a determination of the 
distance with an uncertainty of 20 cm (monitoring also Earth movements at the level of 
centimeters). The real advance of OPERA with respect to previous experiments lies in the 
time tagging system. An upgraded GPS-based timing system at CERN and LNGS allows 
for time tagging with uncertainties at the level of less than 10 nanoseconds. The neutrino 
time of flight is then determined by a statistical comparison between the distribution of the 
neutrino interaction time and the proton probability density function computed from the 
known time structure of the proton beam. The large data sampld^of 16111 neutrino events, 
recorded in a 3-year period, brought the statistical error in the analysis at the same level of 
the estimated systematic error. 

With this procedure, OPERA found the surprising result that neutrinos arrive earlier 
than expected from luminal speed by a time interval pQ 

6t = (60.7 ±6.9 

stat 

±7.4 syst )ns. (1) 

This translates into a superluminal propagation velocity for neutrinos by a relative amount 

6c u = (2.48 ± 0.28 sta t ± 0.30 syst ) x 10" 5 (OPERA), (2) 

where 8c v = (v u — c)/c. The same measurement was previously performed by MINOS (with 
735 km baseline and a broad neutrino energy spectrum peaked around 3 GeV). Although 
not statistically significant, the MINOS result has a central value in the same ballpark of 
the recent OPERA determination [2] 

8c v = (5.1 ± 2.9) x 10" 5 (MINOS). (3) 

In this paper we use the experimental results based on the original data sample reported in the version 1 
of [T] . To remove a possible source of uncertainty, the OPERA collaboration has repeated the analysis with 
a reduced data sample of 15223 events. The final results presented in the version 2 of [T] are consistent with 
those quoted below within one standard deviation. 
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Earlier short-baseline experiments have set upper limits on \Sc u \ at the level of about 4 x 10 5 
in the energy range between 30 and 200 GeV [3]. 

The technical subtleties in the measurement of the neutrino time of flight and the impor- 
tance of this unexpected result call for further experimental scrutiny. Other experiments, 
such as MINOS and T2K, are in a position to repeat the OPERA measurement and this 
can give us decisive information. It is particularly interesting that T2K has a baseline of 
only 295 km. Finding an early neutrino arrival with a time interval of less than half of what 
measured by OPERA would rule out a large class of unknown systematics that may affect 
the present measurement. 

Waiting for further developments on the experimental side, theorists can contribute to 
the debate by analyzing two questions. The first issue is whether the OPERA measurement 
is already in conflict with other tests of special relativity. The second question is whether 
we can envisage simple deformations of the Standard Model that can be reconciled with the 
measurement, possibly suggesting new ways of testing experimentally these scenarios. In 
this paper, we address such questions. 

If the OPERA results are confirmed they will call for a complete reconsideration of the 
basic principles of particle physics. In a Lorentz-invariant theory, propagation of neutrinos 
faster than light would imply violation of causality (i.e. the CNGS beam becomes GS CN 
when seen from a sufficiently boosted reference frame). It is extremely hard to envisage 
a consistent theory having this property. An arguably less radical option is to admit the 
existence of a preferred reference frame that unambiguously defines causal relations but, 
by its very existence, breaks Lorentz invariance. In sect. [2] we summarize and analyze the 
existing constraints on Lorentz violation relevant for our considerations. In sect. [3] we fit 
the OPERA data and other experimental constraints with energy-dependent neutrino limit 
velocities. In sect. [4] we present possible models and their problems. Finally our conclusions 
are contained in sect. [U 

2 Constraints on Lorentz violation 

To assess the compatibility of the OPERA result on 5c u with previous data we should consider 
three sets of constraints on Lorentz violation^] (i) bounds from SN1987a; (ii) bounds on 
flavor- dependent 5c u from neutrino oscillations; (Hi) bounds on Scg, the corresponding limit 
velocity for charged leptons. 

2 For a recent review on experimental constraints on Lorentz violation see ref. [2]. 
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2.1 Supernova SN1987a 



The detection of neutrinos emitted from SN1987a gave us plenty of information not only on 
the process of supernova explosion, but also on neutrino properties. The Irvine-Michigan- 
Brookhaven (1MB) [S], Baksan [G], and Kamiokande II [7] experiments collected 8 + 5 + 11 
neutrino events (presumably mainly u e ) with energies between 7.5 and 39MeV within 12.4 
seconds. The time coincidence of these events can be used to constrain the difference in 
velocity of neutrinos with various energies. Indeed the time delay in the arrival of neutrinos 
with energy E\ and E 2 , respectively, is equal to 

St = ^ [6c u (E 2 ) - SciEj] , (4) 

where g?sn = 51 kpc is the distance of SN1987a. Rescaling a statistical analysis for the case 
of quadratic energy dependence, 8c v oc E 2 [S] , we obtain the bound 

1^(30 MeV) -5^(10 MeV)| < 5 x 1(T 13 at 95% CL. (5) 

However, a larger uncertainty could result from taking into account that the average neutrino 
energy changes with time during the detection interval of about 10 seconds. 

A direct constraint on 5c u , rather than on the difference between 5c u at different energies, 
can be derived from the observation that the optical signal from SN1987a came about 4 hours 
after the neutrino detection, in good agreement with supernova models. This implies 

|<fc(15MeV)| <HT 9 (6) 

This limit should be understood with an order-one uncertainty, which takes into account 
that the precise time delay between light and neutrinos is not known. 

In summary, limits on 5c u from SN1987a are many orders of magnitude more stringent 
than the positive signal of OPERA. One way to evade the SN1987a bound would be to 
assume that Sc u is sizable and positive for neutrinos, but vanishes for antineutrinos. Since 
the detectors that measured the supernova flux were mainly sensitive to antineutrinos, the 
SN1987a observation would not contradict the large 5c u claimed by OPERA. However, we 
are not aware of a concrete setup that realizes this possibility. A more viable solution exploits 
the fact that the SN1987a limits are not directly applicable to the OPERA result, since they 
refer to a different neutrino energy domain, as summarized in fig. [iji. 
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2.2 Neutrino oscillations 



The effective Hamiltonian describing neutrino oscillations can be written as 

2 

771 

H cS = (l + 5c u )E+—. (7) 

Here m 2 is the usual 3x3 neutrino mass matrix and Sc u describes the Lorentz violating 
effect. In general 5c u too could be a matrix in flavor space. The observation that neutrino 
oscillations follow the expected pattern with oscillation phases proportional to L/E, where 
L is the distance travelled by neutrinos, gives strong bounds on flavor-dependent Lorentz 
violations. 

Off-diagonal elements of the matrix 5c ViVj are constrained (9j [TO] and must satisfy the 
bound 

1 n iq / GeV\ /km\ 
M<-»2xl0-(— (8) 

Differences between diagonal elements in the neutrino mass eigenstate basis are similarly 
constrained, because oscillations involving active neutrinos have been observed and the os- 
cillation dip has been clearly seen in KamLAND {y e — > v e disappearance at L ~ 100 km, 
E ~ 10 MeV due to solar u e -> v^ T oscillations [II]) as well as in K2K, MINOS and T2K pj] 
energy spectra [y^ disappearance due to — > v T atmospheric oscillations at L ~ 300 km and 
E ~ 3GeV). Too large neutrino velocity differences would split the wave-packets thereby 
destroying coherence and erasing the observed oscillation dips. Then we obtain the bound 



Sc^-Sc^] <10-^ 21 > , (9) 

where the exact value of the exponent depends on the indices i, j. 

In summary, flavor-dependent Lorentz violations in the neutrino sector are so strongly 
constrained that we can hope to explain the OPERA measurement only with a universal 
limit velocity for all neutrinos. The same conclusion applies to other Lorentz-breaking op- 
erators that affect the neutrino velocity. Therefore, it seems reasonable to assume that 5c u 
is proportional to the identity matrix in flavor space. 



2.3 Charged leptons 

Left-handed neutrinos are part of electroweak doublets. Thus gauge interactions can relate 
Lorentz-violating effects in neutrinos and in charged leptons. This can be a serious concern 
for any attempt to explain the OPERA measurement in a consistent theoretical framework, 
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Figure 1: Top: MINOS, OPERA and FermiLab1979 data on 5c u , together with the upper 
limit from SN1987a (shaded region). Error bars are at la (black) and 2a (blue). Overlaid 
are best fits for 5c u oc E% (dashed red line) and 5c u oc El/(El + E%) (dotted green line). 
Bottom: Upper limits on \Sc e \. The bounds marked with "+" ("—") apply to the case of 
positive (negative) 5c e . See text for an explanation of the different bounds. 
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Figure 2: The Feynman diagram exhibiting how Lorentz violation in neutrinos, represented 
by the red blob, is transmitted to charged leptons £ by quantum electroweak corrections. 

because the present limits on 5c e for electrons are much more stringent than for neutrinos. 
In fig. [TJo we summarize the bounds on 5c e , relating them to the relevant electron energy 
scales. The labels in fig. [T]d correspond to the following observations. 

• "Synchrotron from LEP" is the limit \6c e \ < 5 x 10~ 15 derived from agreement be- 
tween observation and the theoretical expectation of electron synchrotron radiation as 
measured at LEP [13] . 

• "Cosmic rays: e — > e7" refers to the observation that electron vacuum Cerenkov ra- 
diation (the decay process e — > cy) becomes kinematically allowed for E e > m e j \fbc~ e . 
Since cosmic ray electrons have been detected up to 2 TeV, one deduces 5c e < 1CT 13 [T3] . 

• "Cosmic rays: 77 — > e + e~" refers to the process in which high-energy photons are 
absorbed by CMB photons and annihilate into an electron-positron pair. The process 
becomes kinematically possible for E CMB > m 2 e / + 5c e E 7 /2. Since it has been 
observed to occur for photons of about E 7 = 20 TeV, this implies 5c e < 2ml/ E 2 ~ 

10" 15 psj. 

• "Cosmic rays: 7 — > e + e~" refers to the process in which a photon decays into e + e~, 
which becomes kinematically allowed at energies E 1 > m e y/—25c e . Since photons have 
been observed up to 50 TeV one deduces — 5c e < 2 x 10~ 16 [H]. The analogous bound 
for the muon is — <5c M < 10~ n [16J. 

• "Synchrotron from CRAB" refers to the limit — 5c e < 5 x 10~ 15 deduced from astro- 
physical observations of inverse Compton scattering and synchrotron radiation |17j . 

At first sight it may seem that the strong limits on 5c e could be evaded by introducing a 
Lorentz violating interaction with appropriate Higgs insertions able to select only neutrino 
fields. One such example is the space rotationally invariant operator (W L^dolLH), where 
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L is the lepton doublet and SU(2) indices are saturated between the fields in parenthesis. 
Below the weak scale, the Higgs vacuum expectation value leads to an effective operator 
involving neutrinos but not charged leptons. However, quantum corrections due to weak 
interactions generate a non- vanishing 5c e proportional to the original Sc v . This can be seen 
in the previous example, since the operator (H'L*)do(LH) has an anomalous dimension 
mixing with the operator H)L^doL, which contributes to Sc e . 

The effect is quite general and comes from a one-loop diagram with internal exchange of 
a W boson and a neutrino with Lorentz- violating insertion (see fig. [2]), giving 

^"ffv& nw-m (10) 

Barring cancellations, the most optimistic situation occurs when 5c u is as large as the value 
measured by OPERA in a small energy range around -Eopera ~ 20 GeV and vanishes else- 
where. (In sect. [4] we will construct examples of models that fulfill these conditions). In this 
way, from eq. ( JlOj ) we can set a generic lower limit on the value of 5c e based on the OPERA 
measurement: 

5c e > ( E ^^) 2 M^opera) ~ 10~ 9 . (11) 
This result is clearly in conflict with the limits summarized in fig. [T]d. 



To justify eqs. (10) and (11) we consider a simple toy model, with a Yukawa interaction 
between an SU(2) lepton doublet and a sterile neutrino N with Lorentz- violating speed 
c(l + 5cn)- This simplified model illustrates well the dynamics of how 5c e is generated, 
reproducing all the essential features of the more realistic examples later discussed in sect. |4j 
The Lagrangian of the toy model is 

& = £> S m + Ni$N + XHLN + — N 2 (12) 

where k^ = fc M — Sc^u^k ■ u) for any quadri-vector k and = (1, 0, 0, 0). At tree level one 
finds 5c u (k) = (\v/M) 2 5cn up to k < Mj \J5cn, while 5c u is suppressed by 1/k 2 at large k. 

At high energies, much above the W mass, the Yukawa interaction generates at the one- 
loop level a correction to the limit velocity of the left-handed electron (and neutrino) in the 
L doublet: „ 2 

- ■ ., 1 +\n^ + 0(l) 
6(47r) 2 Le A 2 

where we used regularization in d = 4 — 2e dimensions. In the MS renormalization scheme 
the 1/e pole is reabsorbed in the definition of 5cl and the RGE scale p, is identified with the 
low energy scale, Jx ~ Myy. Here A is the maximal scale at which the theory is valid, which 
conservatively could be as low as A ~ M w . 



X 2 5c N 

OC] 



(13) 
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At low energies we can compute the effect by evaluating the contribution of W boson 
exchange using the unitary gauge, such that the only diagram is the one shown in fig. [2] 
Neglecting fermion masses (p 2 = 0), the correction to the electron self-energy S e (p) = fl is 



g 2 (Xv) 2 f d*k lll UhvP L [-g^ + (k + p),(k + p)„/M 2 v ] 



Ml 



(14) 



2 J (2tt) 4 W (k+p) 2 

Here we are interested in the Lorentz- violating term inducing a non- vanishing Sc e , which is 
given by the part of 5E e proportional to i/i. At leading order in Sc^, the neutrino propagator 
acquires a Lorentz-violating term proportional to 



1 



- 5c N 



,k ■ u 



2# 



(k ■ uf 

k A 



and we find a modification of the limit velocity of left-handed electrons, 

^ 2 (Xv) 2 8c N (. 1,1 , u 2 1, 



Sc e (p) 



2(4tt)2 M 2 w , u c ^ w 

The first term inside the curly brackets comes from transverse W polarizations and the second 
(divergent) term comes from the longitudinal W polarization and reproduces the effect in 



1 + 6<7 + h1 ^ 
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1 + Sc n O( 



w 



(15) 



(16) 



eq. (13). The lowest-order term agrees with eq. (11) after taking into account that one needs 



M in the same energy range as -EoperA; while higher-order terms are model-dependent. 
Indeed, in this model we obtain 



Sc e (p < 



w 



Sc e(.P)\uS + 



35 M 2 



;5C U . 



(17) 



VSc^ J vr/IMb 36 (Anv) 2 

This result holds as long asMc Ativ (which is the regime relevant for the models considered 
in sect, [ij because we have used the tree-level relation between 5cn and 5c u . For M 3> Anv, 



the one-loop correction becomes dominant and eq. (17) should be modified 



In the limit of large electron energy (p 2 = with p-u Mw / \J5c^) the effect can be more 
easily extracted from the unbroken SU(2) limit. Since the loop function of the integral in 



eq. (14) depends on the single variable p-u/ 'Mw, the high-energy limit is obtained by setting 
M\y = 0, which corresponds to the calculation in eq. (13). In this regime, for dimensional 



reasons, 5c e remains roughly constant, with a logarithmic dependence on p/Mw, such that 
the effect cannot be fine-tuned away. 



One may believe that the correlation between Sc u and 5c e is specific to the interaction 



in eq. (12). To show that this is not the case, let us consider a different way of breaking 



Lorentz invariance, through a AL = 2 term ^n[Ni{d ■ u)N + h.c], which corresponds to an 
energy-dependent Majorana 'mass-like' term = n(k ■ u). 
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The modification of the neutrino velocity is then given by Sc u = (Xv/M) 2 ^ 2 , where the 
square appears because we consider lepton-number conserving neutrino-to-neutrino propa- 
gation. For the same reason, the correction to the electron limit velocity is also proportional 
to k 2 and consequently again proportional to 5c u , as in eq. (11). Indeed, the only change in 



the loop contribution is that eq. (15) becomes 

1 



where the Lorentz-breaking term has the same form as the last term in eq. (15). Given 



that both Lorentz-breaking terms in eq. (15) give non- vanishing contributions to Sc e with 
comparable order one coefficients and that no other terms are possible, we conclude that the 



estimate of eq. (11) has general validity. 

The muon limit velocity Sc^ is less constrained than in the case of the electron [161 HH] • 
However the strong bounds on Lorentz violation for non-universal flavor interactions prevents 
us from restricting superluminal effects to the v^—y, sector. 

In summary, the existing limits on 5c e present a formidable obstacle for a consistent 
interpretation of the OPERA result. 

Moreover, the unavoidable non-universality due to the charged lepton masses can prop- 
agate at the quantum level an original universal Sc u into flavor- dependent contributions. 



Given the strong bounds discussed in sect. |2.2[ this could also create serious difficulties in 
the constructions of viable theories. 



3 Model-independent analysis 

The OPERA measurement, eq. (|2]), together with the bounds from SN1987a in eqs. ([5])- 
(|6| and those from neutrino oscillations, imply that bc v must be nearly flavor-universal 
and must show a steep energy dependence. It is not difficult to imagine a power law that 
gives 5c u oc E™, with n positive. The problem however lies in the OPERA data themselves 
that show no evidence for an energy dependence in the available energy domain. Using 
only charged- current events inside the detector (a sample of 5489), the collaboration has 
presented an analysis in which the data are split into two bins of comparable statistics, 
containing events of energy larger or smaller than 20 GeV. The average energies of the two 
bins are 13.9 and 42.9 GeV, respectively. The time of neutrino early arrival for the two bins 
is measured to be [lj 

St = (53.1 ± 18.8 stat ) ns for (E v ) = 13.9 GeV 
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Figure 3: Left: Fit of the OPERA and MINOS data assuming 5c„ oc E™. The red contours 
do not include the SN1987a bound (and have Xbest — -U; while the green ones do (with 
x l est = U). The contours are at 1,2 and 3a (68%, 95%, 99.73% CL for 2 dof). Right: 
Same fit assuming Sc u oc E™/(E™ + E™) and including the SN1987a bound. 

5t = (67.1 ± 18.2 stat ) ns for (E u ) = 42.9 GeV. (19) 

The difference between the two values is (14.0 ± 26.2) ns, indicating no significant energy 
dependence. 

The situation is illustrated in fig. [l^,, where the OPERA and MINOS measurements are 
plotted together with the SN1987a constraint. While the large value of 5c u from OPERA 
requires a power index bigger than 2, such a steep energy dependence is in conflict with 
the two-bin analysis. This conflict is quantified in fig. |3^, where we fit the data assuming 
a power law b~c v oc £7™. In the global fit we take into account the energy dependence of the 



OPERA data by appropriately combining the results in eq. (19) 



The red region in fig. [3^, shows the fit obtained including only OPERA and MINOS data: 
the best fit has a very good \ 2 ~ 1 an d a power-law index n fa is preferred. 

The green region shows the global fit, where also the SN1987a data are included. A 
power-law index n < 2 is strongly disfavored. The two regions are statistically incompatible, 
and indeed the global fitj^] including SN1987a has a high \ 2 ~ 14. Furthermore, the global 
fit is incompatible with the older measurements at higher energies [3], shown in figure fig. [I] 

3 Of course only the OPERA collaboration can perform a precise fit. However, although details could 
change, the main features of our analysis are robust. 
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collectively as "FermiLab1979", because of the growth of Sc u at high energy. 

In order to obtain a good global fit, one must assume a distortion from a simple power 
law for 5c v . This situation is less conventional, since it requires the existence of some new 
physics threshold between the energies of SN1987a (about 10 MeV) and OPERA (about 
20 GeV). The same conclusion has also been reached in ref. |19j . 

In sect. [4] we present examples of models that can produce a distortion of the power law, 
as suggested by data. Here, we assume a simple change from power-law behavior at low 
energies to constant behavior at high energies, 5c u oc E™/(E™ + E™). As shown in fig. [3|d, a 
statistically acceptable fit is obtained for n > 3 with E* ~ 10 GeV. MINOS data give a 1-a 
preference towards smaller values of E* and steeper energy growth (n > 4). An example of 
such a fit that well interpolates between data and constraints is shown in fig. where we 
assume n = 3 and 22* = 15 GeV. 



4 Models 

We want to study now how to generate a 5c u with an energy functional dependence different 
than a power law. The safest way is to introduce one or more sterile neutrinos N{ with 
masses Mj within the OPERA energy range and confine the Lorentz-violating interactions 
in this sector. The relevant Lagrangian is 

'Mi 



/A^'WY, ■ + ( '-fNf + XiHLNi + h.c. 



+ -^Lorentz- (20) 



Here J^Lorentz involves only the fields Ni and violates Lorentz invariance. The Yukawa cou- 
plings Aj must be identical for any flavor of the doublet L, but additional contributions to 
neutrino masses may be present in Jz?sm- As sterile neutrinos have no gauge interactions, 
Lorentz breaking will be transmitted to active neutrinos in the lepton doublets L at tree 
level via the active/sterile mixing angles 0j ~ Aji;/Mj <C 1 (experimentally constrained to be 
6>j < 0.1) and to charged leptons at loop level. 

In the low-energy limit (E Mj), the light state mostly coincides with the active one 

flight = v + hNi + 0{E/Mi). (21) 

In order to be compatible with the SN1987A bounds, =£?Lorentz must have zero projection over 
the light state at low energies. Assuming, for example, energy-independent Lorentz breaking 
terms 5cij in the N{—Nj kinetic mixings, one has 

5c v (E < Mi) = J2 OiOj^ + 0(E 2 /M?), (22) 

hi 
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Figure 4: \8cJ\ of the active neutrino (solid red curve) and of the heavy extra sterile neutrinos 
that play no active role (dashed and dotted black curves). The left (right) frame refers to 



the model of section 4-1 (4-%) with the parameter choice described in the text. The dips 
indicate changes of sign of 5c u , where 5c u for the active neutrino is positive for energies in 
the OPERA range and below, and negative at higher energies. 

which vanishes for a sector of two sterile neutrinos where the only non-zero terms are 62 and 
(5ci2- Assuming that only one of the two sterile neutrinos mixes with the active neutrino, 
and that Lorentz-breaking is only present in the sterile neutrino mixing, we automatically 
find that at low energy 5c l/ (E) oc E n with n > 2. The power suppression at low energy can 
account for the absence of anomalous effects in SN1987a. 



4.1 A model with 2 sterile neutrinos 

To protect the SM sector from violation of Lorentz invariance we will adopt the framework of 
ref. (20] based on supersymmetry. In this context no Lorentz violating operators of dimension 
less or equal to four are allowed for the superfields describing SM particles. However, a 
dimension-4 Lorentz breaking kinetic mixing for gauge singlets (such as sterile neutrinos) is 
allowed. In the superfield language this operator has the form 

^Kiju" J d 2 6 Nid^Nj + h.c. (23) 

Here are sterile neutrino chiral superfields with the index i labeling species. The fixed 
vector m m with the components 

^ = (1,0,0,0) (24) 

is introduced for convenience to describe Lorentz breaking in a formally covariant way, where 
space rotational invariance is retained. Finally, the coupling constants form an antisym- 
metric matrix, because the symmetric part of the interaction is a total derivative. Thus, this 



12 



interaction is absent in the case of a single field, and one off-diagonal term is the only possible 
interaction in the case of two fields. Here we consider the case of two sterile neutrinos. 



The operator (23) produces the kinetic mixing for the two sterile neutrinos^ 

^Lorentz = IK W^d^ - IK* tfN^N* . (25) 

Thus the full neutrino Lagrangian including mixings with the left-handed neutrino v is 

if = iva^d^v + iNxa^d^Nx + iN 2 a^N 2 + {^Nl + -^JVf + rmvNx + m 2 vN 2 + h.c.) 
+iKu^N 1 d IM N 2 -iK*u^N 1 d^N 2 , 

(26) 

where m lj2 = Ai j2 t> and we take Mi < M 2 . 

In the high-energy regime (E ^> M\ t2 ), the sterile and active neutrinos decouple and 8c v 
is suppressed. At intermediate energies (Mi <C E <C M 2 ) we can integrate out the heaviest 
of the two sterile neutrinos, finding 

■M l 



if = Wa^u + (mii/JVi + h.c.) + iN l a> 1 d^N l + (-^A 7 ? + h.c.) 



-z/ 2 + i -u^vd^Nx + h.c. 



(27) 



2M 2 M 2 

where we have kept only terms with up to one derivative. There is a non-trivial kinetic 
mixing which, neglecting masses and solving for the equations of motion, gives 

Scu = \S\ 2 , S = ^-. (28) 



M 



2 



The quadratic dependence of 5c u on k arises because the operator of eq. (25) violates lepton 
number, while OPERA observes — > propagation. 

It is interesting to find that the sign of 5c u is fixed and thus neutrinos are necessarily 



superluminal. This is a direct consequence of eq. (25), which predicts sterile neutrinos faster 
than light, while the property of being superluminal is preserved in mixing. The value of 



5c y in eq. (28) is independent of energy, offering the possibility of explaining the flattening 
of the neutrino speed in the OPERA energy range. In practice, however, the intermediate 
energy regime (Mi <C E <C M 2 ) should be very small in order to avoid large mixing angles 
between sterile and active neutrinos. 

In the low-energy regime (E <C Mi,M 2 ) we can integrate out both sterile neutrinos, 
getting the effective Lagrangian up to one derivative 

J? = m»d„v + ( ^^iu^d.p - ^L 2 - + h.c.) , (29) 

M V MiM 2 ^ 2Mi 2M 2 J ' y ' 



4 We work with two-component spinors and follow the conventions of ref . |3T] . The signature of the metric 
is (-,+,+,+). 
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We see that apart from the standard neutrino masses we obtain a Lorentz violating kinetic 
term. However, for a single left-handed neutrino this term is just a total derivative. We 
should then consider higher- derivative Lorentz-violating interactions. Taking for simplicity 
mi = 0, we find 



K 2 m 2 



This leads to a quartic dispersion relation at low energy: 

14 



Thus we obtain 

3|/d 4 |m 2 | 4 £ 2 



E 2 -P 2 - I ^l^ = 0. (31) 



5Cv = 2M 2 M$ ' (32) 



predicting a power-law suppression proportional to E at low energies 



The expression (31) allows us to estimate the energy -E tra ns where the behavior of 5c u 
changes from a quadratic growth to a constant. Equating 



1*1 



i 



E ~ \6\ 2 E 2 , (33) 



Ml 
we find 

^trans ~ M X /\5\ . (34) 

This shows that the transition occurs at energies parametrically larger than M\. 

These analytic estimates are confirmed by a numerical calculation in which we compute 
exactly the dispersion relation Ei = Ei(p) for the three neutrino states derived from their 
Dirac equations. The neutrino group velocities is obtained by taking the derivative of Ei 
with respect to momentum: Cj = dEi/dp. The numerical result corresponding to the choice 
M 2 = 2Mi = 40 GeV, 6 2 = 0.1, k = 0.7 is shown in fig. |§l 

This model has a potential problem. The contribution to active neutrino masses in 



eq. (29) obtained by integrating out N% and N 2 is too large. Of course this can be cancelled 
by other contributions to Majorana masses of active neutrinos. Although this may seem 
to be a special fine tuning, the cancellation can occur naturally. One example is the case 
in which the active neutrinos have Yukawa couplings to a single sterile heavy neutrino N$, 
while N 12 mix with N 3 in the mass matrix. This generates a rank-one mass matrix for the 
lighter states. So this problem can be easily bypassed. 

Another difficulty is that the model produces at low energy (E <C Mi^) an effect sup- 
pressed only by E 2 , which is not enough to fit OPERA compatibly with the SN1987A 
bounds, as illustrated in fig. [4^,. This problem can be solved in different ways, as shown by 
the following examples. 
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4.2 A model with quadratic mixing term 



To remedy to the E 2 low-energy behavior of 5c u we can add a further derivative in the 



Lorentz-violating term in eq. (|25|) and use the Lagrangian 

(M 
\~2 



= ivo^d^v + iN^d^Nx + iN 2 a^d^N 2 + 
+Ku u u»d u N 1 d l ,N 2 + k* vfu^dvNxd^Nz . 



fMl Nl + ^Nl + m 2 uN 2 + h.c. 



(35) 



In the high-energy region (E 3> M 2 ) 5c u is suppressed and dies out. In the intermediate 
region Mi <C E <C M 2 (or more precisely i?trans "C E <C M 2 , where -Etrans will be computed 
in eq. (42) below), we obtain the effective Lagrangian 

Km 2 



Lorentz 



M 2 



«V vdyduNi + h.c. 



(36) 



where we have kept terms up to two derivatives. This produces the dispersion relation 



E 2 -p 2 - 



\k\ \m 2 \ 



-£ 4 = 



Hence 



3\K\ 2 \m 2 \ 2 E 2 
2M? 



8c v 



(37) 



(3* 



We see that in this case the velocity of neutrinos is not constant in the intermediate range 
but depends quadratically on energy. 

At low energies (E <C Mi) we integrate out the remaining sterile neutrino and obtain 



This yields the dispersion relation 



K 2 m 2 r, 



2M 1 M 2 ^ " ^ 



E 2 -p 2 



\K \ \m 2 \ 



M{Mi 

and hence a Sc u scaling as the sixth power of energy. 



E s = 



8c v 



7\K\ A \m 2 \ A E G 
2M 1 2 M 2 4 



(39) 



(40) 



(41) 



By comparing eqs. (40) and (37) we obtain that the transition between the two regimes 
occurs at 

/ M X M 2 , . 

(42) 



EL 



nm 2 



These analytic estimates are confirmed by our numerical calculation, as exemplified in 
fig.|3>, which corresponds to the choice M 2 = 1 TeV, Mi = 1 GeV, 6 2 = 0.02, « = (30 GeV)" 1 
(in the two-derivative case k is a dimensionful coupling). In this case we obtain a satisfactory 
fit of data and constraints, because of the steeper dependence on E at low energies. 
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4.3 A model with three sterile neutrinos 



An alternative approach is to employ three species of sterile neutrinos with Lagrangian 

Mi 

i<j 

where i,j = 1, 2, 3. We assume the hierarchy of masses Mi <C M 2 <C M 3 . 



Se = m^u + m^Ni + [rmvNi + ^N? + ]T /K,,//".V,r7 / ,.V / + h.c.) , (43) 



To simplify the calculations let us assume mi = m 2 = «i2 = 0. At large energies 
E 3> M 3 the active neutrino decouples from the steriles and propagate with the speed of 
light. At M 2 <C E <C M 3 we integrate out the third neutrino and obtain the following 
Lorentz-breaking contributions 

O, K 13™3 • a o »r , K 23™3 . „ o Ar 



» 1 M 3 



(44) 



The velocity of the active neutrino in this regime is given by 

X |jf |2 , Ijf |2 X _ K 13 m 3 jj _ ^23^3 ,,_v 

5c u = \5 13 \ + \5 2 3\ , di 3 = vir , O23 = vir • (45) 

M 3 M3 

At lower energies (M\ <C E <C M 2 ) we integrate out the second neutrino and the Lorentz 
breaking part becomes 

-2 ,J2 ™2 



where we have considered only terms up to two derivatives. The third term here gives the 
quadratic growth of Sc u with energy. 



Finally, at E <C Mi we integrate out N\ and obtain 



y> m 



2M% 



K 13 _|_ ^23 

Mi M 2 



K<V) 2 . (47) 



The quadratic growth can be shut off at the price of the fine-tuning condition (ki 3 /k 23 ) 2 — 
— Mi/M 2 . The next operator is quartic in derivatives and will give an E 6 dependence of 8c v 
on energy. 

Let us conclude our discussion on these models with an important consideration. The 
models described here realize the optimal situation in which 5c u is large only at OPERA 
energies: at lower energies there is a power-law suppression and at higher energies bc v de- 
creases because energy dominates over the active/sterile mass mixing term. In spite of that, 



16 



the quantum corrections to the velocity of charged leptons £ is too large. For instance, in 



the model of sect. 4.1, we find 



See 



A2 



(4tt) 



(4f 



in agreement with the generic estimate of eq. (11). Note that supersymmetry cannot be 
invoked as a solution, because at momenta of the order of the OPERA energies the wino 
loop contribution cannot cancel the W ± loop. 



5 Conclusions 

It has been often speculated that Lorentz invariance can be violated in new frameworks 
beyond the Standard Model, especially in the context of various quantum gravity scenar- 
ios, such as space-time foam [22], string theory [23], Hofava gravity [21], extra dimensions 
with asymmetric space-time warping factors [25], or inhomogeneous Higgs backgrounds [26J. 
Tests of Lorentz invariance can provide us with precious hints of physics at energy scales 
much beyond those directly probed in colliders and the experimental discovery of an actual 
Lorentz violation in nature would have a revolutionary impact on our understanding of the 
fundamental principles of the particle world. 

The recent OPERA measurement, pointing towards evidence for superluminal propaga- 
tion of neutrinos, opens an exciting perspective. Careful scrutiny of the OPERA methodol- 
ogy and confirmation from other experiments is certainly required before the measurement 
can be accepted. From the theoretical point of view, the result is quite surprising, not be- 
cause of the violation of a sacred principle of special relativity, but because of the size of the 
effect. The amount of Lorentz violation in the neutrino sector required to explain the effect 
is much larger than the typical bounds that have been extracted for other Lorentz-violating 
parameters. This poses a question of logical consistency of the new experimental result. 

In this paper, we have shown that the OPERA measurement is not necessarily inconsis- 
tent with other tests of special relativity in the neutrino sector, but requires very particular 
assumptions on 5c u , the relative difference between the neutrino and photon limit velocities. 
Various observational constraints, reviewed in sect. [2j force 5c u to be nearly flavor universal, 
to drop steeply at low energies (E v < 10 GeV) and remain roughly constant in the OPERA 
energy range (around 10 to 40 GeV). Such a peculiar behavior is unexpected but, as shown 
in sect. [4] with the help of several illustrative models, can be accounted for by introducing 
sterile neutrinos with Lorentz-violating interactions and masses in the OPERA energy range. 

As discussed in sect. [2j a fairly generic and formidable constraint on any interpretation 
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of the OPERA measurement comes from Sc e , the relative difference between electron and 
photon limit velocities. Even if we restrict the original Lorentz violation to the neutrino 
sector, quantum corrections due to weak interactions will generally pollute the electron 
sector as well. Barring unknown cancellations, the resulting effect is too large, showing 
conflict between present bounds on 5c e and the OPERA measurement. Fine-tuning the 
value of 5c e at a particular energy cannot cure the problem, because there are strong bounds 
on 5c e applicable at several different energies. 

Much work has still to be done, both on the experimental and the theoretical front, before 
the fate of the OPERA measurement and its impact on physics will be determined. 



Note added 

After this paper was published, the authors of ref. [27] presented another strong constraint 
on superluminal interpretations of the OPERA result, based on energy loss from neutrino 
decay. Their constraint is as numerically stringent as the constraint presented here, based 
on loop contributions to 5c e , but more direct. Ref. [2Zj assumed an energy-independent 5c, 
but their argument is fairly robust, since any dispersion relation leading to a superluminal 
neutrino, like the ones we consider, implies that a neutrino decay process is kinematically 
allowed. Indeed, let us consider the decay process u{p) — > v(p')X, where X = e + e~ or any 
light final state. The process is allowed when the invariant mass of X is positive or, in other 
words, when the difference of the 4-momenta of the two neutrinos is time-like, 

E(p) - E{p) >p-p. (49) 

Here E(p) is the neutrino energy for a generic dispersion relation and p (p') is the momentum 
of the initial-state (final-state) neutrino. For a superluminal neutrino we approximately find 

<LE E(p) - E(p') , . 

C >1 =4> c= — — > 1. 50 

dp p — p 

Thus a superluminal neutrino always has kinematically-allowed decay channels, indepen- 
dently of its dispersion relation. 
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